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Abstract 

The effect of Cooperons on metal-insulator transitions (MIT) in disordered in- 
teracting electronic systems is studied. It is argued that previous results which 
concluded that Cooperons are qualitatively unimportant near the MIT might 
not be correct, and that the problem is much more complicated than had 
previously been realized. Although we do not completely solve the Cooperon 
problem, we propose a new approach that is at least internally consistent. 
Within this approach we find that in all universality classes where Cooperons 
are present, i.e. in the absence of magnetic impurities and magnetic fields, 
there are logarithmic corrections to scaling at the MIT. This result is used 
for a possible resolution of the so-called exponent puzzle for the conductivity 
near the MIT. A discussion of the relationship between theory and experiment 
is given. We also make a number of predictions concerning crossover effects 
which occur when a magnetic field is applied to a system with Cooperons. 
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I. INTRODUCTION 



The current theoretical description of the metal-insulator transition (MIT) problem as- 
serts that the presence or absence of the particle-particle or Cooper channel does not quali- 
tatively modify the MIT.S'i At first sight this may seem surprising, as much of the early work 
in the modern (post- 1979) theory of the localization problem concentrated on the Cooper 
channel and the backscattering or weak localization effects it produces.! Also, numerous 
experiments confirmed the presence of these weak localization effects in weakly disordered 
metallic systems in the absence of magnetic impurities and magnetic fields.§ However, the 
assertion appears less surprising if one recalls that electron-electron interaction effects in 
the presence of disorder lead to many of the same effects as Cooperons.! If this is so in 
the weak disorder regime, and if one acknowledges that electron-electron interactions are 
in general relevant for the MIT, then it is conceivable that Cooperons do not lead to any 
additional effects at the MIT over and above those produced by the interplay of interactions 
and disorder alone. This was in fact the conclusion reached by Finkel'steinlll and other^ 
who have argued that the Cooper channel is irrelevant, in the sense of the renormalization 
group (RG), for the MIT, that interaction effects effectively replace Cooperon effects near 
the MIT, and that MIT with or without Cooperons are qualitatively the same. 

In this paper we argue that the latter conclusion is probably not correct. We present 
a consistent description within which the Cooper propagator, or the effective Cooper inter- 
action amplitude, r^, is a marginal operator rather than an irrelevant one. This marginal 
operator leads to logarithmic corrections to scaling that are characteristic for those univer- 
sality classes where Cooperons are present. 

Technically, we first show that the MIT problem in the presence of Cooperons is substan- 
tially more complicated than had previously been realized. In particular, the renormalization 
procedures used in existing treatments of the problem do not lead to a finite renormalized 
field theory. We further argue that even if this problem is ignored, the RG fixed point 
structure obtained in previous treatments is not stable with respect to the consideration of 
higher order terms that were neglected. We then present a partial solution to this problem. 
We first point out that a principal unanswered question is that of how many renormalization 
constant are needed to make the theory finite. This is closely related to the physical ques- 
tion of whether Tc is a simple scaling variable near the MIT or whether it consists of several 
scaling parts. Since at present we do not have a firm answer to this question, we first derive 
RG flow equations for all other coupling constants which appear in the theory. These flow 
equations can be expressed in terms of whose scaling behavior is a priori unknown. We 
then derive an Eliashberg-type integral equation for F^. General arguments lead us to the 
conclusion that F^ consists of several scaling parts, and that it approaches its asymptotic 
value at the MIT, F*, logarithmically slowly. 

The most important conclusion from our considerations is that in systems without mag- 
netic impurities or magnetic fields there are logarithmic corrections to scaling at the MIT. 
This in turn implies that it is virtually impossible to experimentally reach the asymptotic 
critical scaling regime for these universality classes, and that the existing experiments mea- 
sure effective exponents, rather than asymptotic critical exponents. We show that these 
logarithmic corrections to scaling provide a possible resolution of a long-standing problem. 
The critical exponent s for the conductivity in Si:Pi and some other systems,0il all of which 
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are believed to be in universality classes that contain Cooperons, is experimentally observed 
to be smaller than 2/3, in apparent violation of a rigorous bound that requires s > 2/3 in 
three-dimensional (3-D) systems.! We will see that logarithmic corrections to scaling can 
easily account for an apparent or effective exponent Sg// ~ 0.5 even if the true asymptotic 
critical exponent obeys s > 2/3. The same conclusion was reached in a previous short 
account of part of the work presented in this paper .0 

The plan of this paper is as follows. In Sec. ^ we use Finkel'stein's effective field theory 
for the to calculate the perturbative corrections to the coupling constants which appear 
in the theory to one-loop order. In Sec. |T| we use a normalization point RG approach to 
renormalize the field theory and to derive RG flow equations for the coupling constants in 
terms of the Cooper propagator Tc. We also show how previous results can be obtained 
from certain assumptions and approximations for Tc- We find that these assumptions do 
not lead to a finite renormalized theory, and that terms neglected in the previous treatments 
invalidate this approach in any case. In Sec. |IV| we classify and discuss the possible solutions 
of the integral equation for Tc- We argue that generic solutions of the integral equation yield 
a Tc that approaches its fixed point value logarithmically slowly, which leads to logarithmic 
corrections to scaling. In Sec. we discuss the experimental consequences of this result. We 
compare the theory with existing experiments and suggest a number of new measurements 
to further test the theory. In particular we discuss crossover effects due to an external 
magnetic field which changes the universality class and eliminates the logarithmic corrections 
to scaling. We conclude in Sec. with a general discussion of our results as well as the 
current status of the MIT problem. 

II. THE FIELD THEORY, AND THE LOOP EXPANSION 

In the first part of this section we recall the basic field theoretic description of the 
disordered electron problemSfl and derive the Gaussian propagators of the field theory. We 
then explain how to obtain perturbative corrections to the coupling constants by considering 
the generating functional to one-loop order. 

A. The Model 

The existing theoretical description of the MIT is based on the usual assumption of the 
theory of continuous phase transitions, viz. that the physics near the transition is dominated 
by the low-lying excitations of the system. As for the description of other phase transitions, 
the problem is then reduced to obtaining a solution for an effective field theory for these slow 
modes. Once the field theory has been identified, the technical apparatus used to obtain a 
solution is the RG. The initial problem is to obtain the field theory, which requires a physical 
identification of the relevant slow modes. 

In the field theory which we will use,0'i the slow modes in the metallic phase are assumed 
to be the diffusion of mass, spin, and energy density. The effective field theory for long- 
wavelength and low-frequency excitations that describe how these diffusive processes change 
across an MIT is defined by an action. 
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S[Q] = -^Jd^ tr(vg(x))' + 2H Jd^ ^r(^]g(x) 

-— E ^«yrfx[Q(x)7(")Q(x)]. (2.1) 



Here the field variable Q is an infinite matrix whose matrix elements, Q"^, are complex 4-by- 
4 matrices (spin-quaternions) which comprise the spin and particle-hole degrees of freedom. 
The labels a, P = 1,2, - ■ ■ ,N denote replica labels. In deriving Eq. (p.l]), quenched disorder 
has been integrated out by means of the replica trick,lli and the limit — > cxo is implied at 
the end of all calculations. n,m = — oo, ■ ■ ■ , +oo are Matsubara frequency labels. Q = lun 
with I the identity matrix, ujn = 2TTT{n + 1/2), is a fermionic frequency matrix, and tr 
denotes a trace over all discrete degrees of freedom. [Q'^^'^'^Q] in Eq. ( |2.1| ) is defined as. 



W^Q]= E 5„,+„3,n2+n.E E (-!)'■ 
nin2n3?i4 a r=0,3 

xtr((r,®.o)Qrn2)^K(^^'®"o)^"""J ' (2-2a) 



[Ql^'^Q] = - E ^n,+n3,n2+n, E E ("1)^ 
nin2n3ri4 « r=0,3 

xj2tr{{rr ® s,)QT,^^)tr{{Tr ® s.)Q°^„J, (2.2b) 

i=l 



EE 

n\n2nzn4 a r=l,2 

tr.((r. ® So)Q^^rn2)^^-((^'- ® ^o)Q::^„J] , (2.2c) 



with tr = trgtrj. where tr^- acts only on the r's and iVg acts only on the s's. In Eqs. ( p. 21 ), 
'i'o = So = o"o and = — = —icTj {j = 1,2,3), with aj the Pauli matrices. The theory 
contains five coupling constants. G = S/naB with as the bare or self-consistent Born 
conductivity is a measure of the disorder, and H = TrNp/A plays the role of a frequency 
coupling parameter with Np the bare density of states (DOS) at the Fermi level. Kg and 
Kt are singlet and triplet particle-hole interaction constants, respectively, and is the 
singlet particle-particle or Cooper channel interaction constant. At zero frequency the triplet 
coupling constant in the particle-particle channel vanishes due to the Pauli principle. A 
disorder generated, frequency dependent triplet particle-particle interaction constant has 
been discussed elsewhere. El For simplicity we formulate the theory with a short-range model 
interaction, i.e. the Kg^t c are simply numbers. For the more realistic case of a Coulomb 
interaction Kg is x-dependent and must be kept under the integral in Eq. (|2.1|). Most results 
for this case are easily obtained after all calculations have been performed by essentially 
putting K. = -H,B but occasionally subtle complicafons arise. Since these are purely 
technical in nature, and decoupled from the Cooper channel induced problems which are 
the subject of this paper we will not discuss them. We will, however, give results for the 



Coulomb interaction case in Sec. Ill 



The matrix Q is subject to the nonlinear constraints. 
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(2.3a) 



trQ = , (2.3b) 

and satisfies, 

Qt = C^Q^c = Q , (2.3c) 

with, 

C = t{n ® S2) . (2.3d) 

The Ti are the quaternion basis and span the particle-hole and particle-particle space, while 
the Si serve as our basis in spin space. For convenience we expand Qnm ^^is basis, 

Q:i = j:j:iQniirr^s.) . (2.4) 

r=0 i=0 

From the first equality in Eq. ( p.3c| ) it follows that the elements of Q describing the particle- 
hole degrees of freedom are real, while those describing the particle-particle degrees of free- 
dom are purely imaginary, 

rQnm ~ rQnrn 1 (^ ~ 0; 3) , (2.5a) 

rQnrn rQnrn 5 ('^ ^i"^) 5 (2.5b) 

In addition, from the hermiticity condition one obtains, 

°Ql = (-ir°Q^n , (r = 0,3) , (2.6a) 



QZ = i-iy-''lQZ , (r = 0,3; z = 1,2,3) , (2.6b) 

rQfm = rQi"r^ , (^ = 1, 2) , (2.6c) 



IQfm = -IQ^rl , (r = 1,2; ^ = 1,2,3) . (2.6d) 

The constraints given by Eqs. ( |2.3a) , F^^^ and the hermiticity requirement, Eq. ( |2.3c| ), can 
be eliminated by parametrizing the matrix Q by,0 

(1 - ggt)i/2 _i for n > 0, m > 

^ _ I q for n > 0, m < (o 7 ^ 

gt forn<0, m>0 ' ^ 

-(1 - gl'g)^/^ - 1 for < 0, m<0 



Here the q are matrices with spin-quaternion valued elements g"^; n = 0, 1, 
— 1, —2, ■ ■ -. Like the matrix Q, they can be expanded as 



m 



5 



<l:^=j:j:l€^irr®s.) . (2.7b) 

r=0 i=0 



Note that the q do not satisfy the symmetry relations given by Eqs. (p.6| ) for the Q. 

We have so far given the theory for the so-called generic (G) universality class which is 
realized by systems without magnetic fields, magnetic impurities, or spin-orbit scattering. 
Apart from class G, the second universality class with Cooperons is the one with strong 
spin-orbit scattering (class SO). For class SO the action is shown as above, except that the 
particle-hole spin triplet channel is absent, i.e. the sum over the spin index i in Eqs. ( [2.4| ) 
and ( |2.7b|) is restricted to i = O.El In what follows we will give results for both class G and 
class SO. 

With the help of Eqs. (|2.7| ) one can expand the action in powers of q. 



S[Q] = T.Sn[q] , (2.8) 

n=2 

where Sn[q] ~ q"'- We first concentrate on the Gaussian part of the actionjll 

S2[q] = -4 / E E ;?i2(p);Mi2,34;g34(-p) . (2.9a) 

•^P r,i 1,2,3,4 

Here /p = / (ip/(27r)'^, and 1 = {ni,ai), etc. The matrix M is given by, 

^,3Mi2,34(p) = ^i^^{5i3[p' + GHiuj,,, - un,)] + 5 ^,^,5 G K , (2.9b) 
where = s, 1^1^2,3 = t, and 

i,2^i2,34(p) = -^^^{SisIp^ + GH{ujn, - uJn,)\ + 5„,„3 (5,,o27rTGir J , (2.9c) 
The Gaussian propagators are given by, 

^gi2(pi)ig34(p2))^'^ = / D\q\ ;gi2(pi) ig34(p2) e^^iy I D\q\ e^^M 

= ^5,,5,,(27r)^5(pi + P2);Mf2;34(Pi) , (2.10a) 



with, 



0,3Mr2;34(p) = 5l_2,3-4G(5i3P„,-„,(p) + ^V^^^_^^{j>) ) , (2.10b) 



ni-n2 

{p)-6 



G2'kTK^ 



l + G2'KTKJn,+nM 

where. 



(2.10c) 
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Up) = E (P) ■ (2.10d) 

ni>0,n2<0 

Here we have introduced the propagators, 

Vn{p) = [p^ + GHnn]-' , (2.10e) 

= [P' + G{H + Ks,t)nn]-' , (2.10f) 

AD:'*(p)=P:-*(p)-P„(p) , (2.10g) 

with Qn = 27[Tn a bosonic Matsubara frequency. Physically, and are the energy, 

mass, and spin diffusion propagators.0 Equations ( ^■10b| , p!lOc| ) can be put into a more 
standard form by summing over, e.g., and 71,4, 

E 0,3^r2!34(p)/G'=(l-5«,«2)^^ni-n2(p)+'5aia2^^n;-n2(P) > (2-lla) 
n3>0,n4<0 

E l,2^r2!34(p)/G' = -(1 - 5aia2)^^n,i-n2(p)-^aia2(l " '^io)^^ni-n2 (p) 
n3>0,n4<0 

-^»0^aia2 . , ^ . (2.11b) 

1 + G2'nTKcfn^+n2{p) 



Examining the various terms in Eqs. ( |2.11| ) we see that all of them have a standard propaga- 



tor structure except for the last contribution in Eq. (p.llbj) . In interpreting these propagators 
as having a standard structure, the Matsubara frequencies in Eqs. ( |2.10e| )-( |2.10gD are taken 
to be analogous to a magnetic field at a magnetic phase transition, i.e., the MIT occurs at 
fi„ — and Vtn or the temperature is a relevant perturbation in the RG sense. Using Eq. 
( p.lOap , changing the sum to an integral, and placing an ultraviolet cutoff, f2o, on the re- 
sulting frequency integrals shows that 27rT/„j+„2(p) in Eq. (|2.11b| ) diverges logarithmically 
in the long wavelength, low temperature limit. With i^'c > 0, we see that the last term in 
Eq. ( |2.11b|) is logarithmically small compared to the other terms in Eqs. (|2.11| ). 

We conclude this subsection with two remarks. Firstly, the logarithm discussed above 
that appears in the particle-particle density correlation function is just the usual BCS log- 
arithm. However, since we consider a system with a repulsive Cooper channel interaction, 
Kc> ^1 which is not superconducting in the clean limit, this does not lead to a Cooper insta- 
bility. Rather, the last term in Eq. (|2.11b|) vanishes logarithmically in the limit p, T — 0. If 



the structure of this term persists for disorder values up to the MIT, and if it couples to the 
physical quantities like, e.g., the conductivity, then it will lead to logarithmic corrections to 
scaling. Secondly, considering the Gaussian theory one can already anticipate a fundamental 
problem with any RG treatment of the field theory. To see this, note that at the Gaussian 
level the two-point vertex functions are given by 6*2 [g]. Examining the Eqs. ( |2.9| ) we see that 
at this order no singularities, neither in the ultraviolet nor in the infrared, are present in the 
vertex functions, and there is no explicit cutoff dependence. This should be contrasted with 
the corresponding two-point Gaussian propagators given by Eqs. (|2.10| ). Because of the last 



term in Eq. (|2.10c|) , both an infrared singularity and a dependence on an ultraviolet cutoff 
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appear. In the usual RG approach such cutoff dependences are ehminated from the field 
theory by the introduction of suitable renormalization constants. Here, unusual features are 
that vertex functions and propagators behave differently with respect to their cutoff depen- 
dence, and that the cutoff dependent term is logarithmically small rather than large. In 
previous RG treatments of this problem, the procedure used was effectively to renormalize 
Kc in Eqs. (|2.9c| ) and ( p.lOc|) such that the renormalized two-point propagators were finite 
as f2o — cxo. It is easy to see that such a procedure leads to renormalized two-point vertex 
functions that contain a singularity in this limit, cf. Eq. ( |3.10|) below. In Sec. |Tl| we will 
further discuss this problem, and we will propose an alternative renormalization procedure. 



B. One-loop perturbation theory 



We will not be able to present a final solution of the Cooperon renormalization problem, 
we will therefore consider several renormalization procedures, and will discuss which ones 
are at least internally consistent and which ones are not. One of these procedures is based 
on perturbation theory for the generating functional for the vertex functions, r[g], which is 
closely related to the thermodynamic potential. The perturbation expansion for r[g] can be 
generated by standard techniques.^ First, a source term with an external potential J"/^(x) 
is added to the action so that the average of g"^, which we denote by g"^, is nonzero. We 
can then obtain a loop or disorder expansion for T[q\. It is convenient to expand T[q\ in 
powers of g. 



°° 1 r 

= I] T7y d^i---d^N r^3,4.....2iV-l,2Ar(xi, " " " , X^v) gi2(xi) ■ ■ ■ g2JV-l,2Jv(x7v) 

N=l ■ l,2,...,2Af 



where F*^^) is the N-point vertex function. 
To zero-loop order, one has 



Tofel = -S[q] 



(2.12) 



(2.13) 



At higher-loop order the various coefficient in S acquire perturbative corrections, and in 
addition structurally new terms are generated. For simplicity we restrict our considerations 
to the two-point vertex function F^^-*, and to the one-point vertex function F*^^-* which is 
related to the one-point propagator P^^^ =<o Qnni^) >• one-loop order F^^^ is given by 
the second derivative with respect to q of the right-hand side of Eq. ( p.9a| ) with q ^ q, and 
with the matrix M replaced by 



0,3^(2,34 (P) 



Si- 



■2,3-4 



and. 



13 



p 



G 



UJ. 



n2) 



n\n2 



(2.14a) 



1+2,3+4*1 ^13 



p 



Gnir 



n4 



(2.14b) 
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where Gnma, -^nma, and K^^'^^^^^^^ are given by G, H, and Ks,t,c, respectively, plus frequency 
dependent one-loop perturbative corrections. In our notation we have suppressed the fact 
that these corrections are in general also momentum dependent. In the absence of Cooperons 
these corrections have been discussed in detail elsewhere.^ In general the momentum and 
frequency dependence of the corrections is quite complicated. Here we just give the results 
to leading order in 1/e, with e = D — 2. With A an ultraviolet momentum cutoff, and 
G = GSD/i^ir)^ with So the surface of the D-dimensional unit sphere, we obtain. 



G 

Gn\n2 G -\- G — 

4e 



'n2/ 



G^ / [/i^(p,l]„,,_„J + 3/l(p,l]„,_ 



+/^(p, + G^ f [/2^(p, n^.-n^) + 3/^(p, f^„,-„J] , (2.15a) 

Hn,n2 = H + GH |^[/,^(p,fi„,_„J+3/l(p,fi„,_„J + ^/r(p,f^„J + ^/^(p,^]-„J] 

-GH [ [J^ip, + 3 J*(p, , (2.15b) 
Jp 



K. 



n2 } 



(2.15c) 



K 



t 



-G2K, I /^(p,fi„,_„J - ^{H + K^f I J2^(p,fi„,_„J , (2.15d) 



^n.n^nan, = - ^[A^ " (Gi/r]„3_„ J^/^] (TT, - SK^) +G [ r,{p, , (2.15e) 



Here, Qq = 0{ls? /GH) is a cutoff frequency, and. 



m 



(2.16a) 



/r(p,fi. 



--G27rT y ^^^^ r^f^^m(P 

8 1 + G27rTir J_„(p) V 



(2.16b) 



/-(2 no/27rT 

/|'*(p,fi„) = - — (ir,,,)227rT E fi^P:;*(p)(p„+^(p 



X 



1 - 2p2l)„+™(p) , (2.16c) 



1 Ott'V ^o/27rT 

jf(p,fi„) = -(K,,/i/)-— E [AP^:;*(p)+GKm^^^^^*(p)P^+„(p)] , (2.16d) 



n m=n 
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^7cT ^ ^ r TTl Tfl ~\ 

Ji(p, = ^ (1 - -)P^+„(p) + -P^(p) , (2.16e) 

m=l ^ ^ 



f2o/27rT Q,Q/2-nT ^ Q.a/2-KT 

J3(p,^^n) = - E -APj„(p) - GiTi E I^m(p)2^n^+n(p) + -i^. E ^^^^ (p)^^n^+n (p) 



+T^* E ^^^:.(p)^^„.+n(p) + -i^* E (2-^)i^J,Xp)A.+.(p) , (2.16f) 



^ ^ 1 + AcGivri /_m(,pj 



1 _ 1 . _ . ii'. 



Jr(p,fi„) = -G27rT E -AP:,(p)— -— -^r^— - , (2.16h) 
For the one-point vertex function one finds, 

^«(^]„) = l + G' / [/r(p,i^„,) + 3/J(p,fi„) + Ji^(p,fi„,)] . (2.17) 

■J D 



In giving Eqs. ( |2.15| )-( |2.17| ) we have neglected terms that are finite in D = 2 as (A, fio) ^ oo, 



and a delta function constraint (cf. Eqs. (|2.14|) ) is understood in Eqs. ( p.l5c| , p~15d|J2.15e]) . 



Some of these terms depend on ni and n2 separately, not just on the difference ni — n2. In 
Eq. (|2.15b| ) we have written a separate dependence on ni and n2 explicitly for later reference. 
Also, the complete frequency dependence of K^^-^^^ nam ^nin2 nam more complicated 
than the one shown. However, for most of our purposes it is sufficient to treat all 'external' 
frequencies as equal, and we do not have to deal with the (substantial) complications that 
arise from the full perturbation theory. 

Finally, let us discuss one important point. To one- loop order the two-point propagators 
are given by the right-hand side of Eq. ( p.lOaj ) with the inverse matrix replaced by the 



matrix M'~^ which contains the perturbative corrections. For the particle-hole degrees of 
freedom one can show in general that, except for irrelevant terms, the matrix M'~^ has the 
same form as the matrix M~^, with the only difference being that the corrected coupling 
constants appear in M'~^. The underlying reason for this feature is the conservation laws 
for mass, spin, and energy density. For the particle-particle degrees of freedom the situation 
is different. There are no conservation laws which guarantee that the form of the last term 
in Eq. ( p.lOc| ) will not change at higher orders in the loop expansion. In general, the matrix 



i_2(M') is given by Eq. (|2.10c|) with the replacement. 



l + G27rT£/„,^_„,(p) ^^--'"^-^^^ ' ^^-^^^^ 
where F'^'^^ satisfies an Eliashberg-type integral equation, 

nin2,n3n4 nin2,nLn' ^2 In ,, i W ,,(',, , i , , \ ^^nin2,n3n4, ■ V^--^"^// 
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Here we have again suppressed the momentum dependence of V^. Note that if we make the 
substitutions K^,^, ^^^^ — >• i^c, -f^n'^n'j ~^ H, Gn^n'^ — >■ G, replace the remaining sum by an 
integral, let p, T — >• 0, and use an ultraviolet cutoff Qq on the frequency integral, then F'^ 
has the standard BCS form, 

r(p,w„,+„J = — , (2.18c) 

with 7^°) = KJH, and D'-^^ = 1/GH. 

The actual Cooper propagator to all orders in perturbation theory would have the simple 
form given by Eq. ( 2.18c| ) only if the coupling constants K'^ and H were constant to all orders. 
In general, the structure of the Cooper propagator will be more complicated and to obtain 
it one has to solve the integral equation, Eq. ( |2.18bD . For later reference we symmetrize the 
equation by defining 



[-'^ n 1 712 na n4 J 

and 



F'^ 

F — nin2,n3n4 /r, i qi N 



3n4j 



For T ^ the integral equation for F then reads. 



i [uj,\l,uj ) + I duj - 



D{uj',-n -uj')p^ + 2uj' + n 

= -f{iu,n,uj") , (2.20a) 

with, 

D{uj', -n - Lj') = [G{iu', -n - uj')H{uj\ ~n - lj')]'^ , (2.20b) 

Since we consider the zero temperature limit we have made the replacements, 

27rT(ni + = 2-nT{n^ + n^) —VL 

2TiTn\ uj' 
2'KTni —>■ uj" 

f „i„2,„3„4 T{uj, ^, uj") (2.20c) 

etc., and for definiteness we have assumed i7 > 0. We stress again that Eq. ( p.20a|) has the 
form of an Eliashberg equation with a repulsive kernel. 



III. THE RENORMALIZATION GROUP FLOW EQUATIONS 

In the first part of this section we review the general procedure of the field theoretic RG 
as applied to a nonlinear sigma-model. We then use a normalization point RG procedure to 
obtain RG flow equations for all of the parameters that appear in the field theory except for 
the Cooper propagator F. Finally, we discuss the reasons why previous attempts to derive 
a fiow equation for F are problematic. 
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A. The Nonlinear Sigma-Model and the RenormaHzation Group 

The basic philosophy behind the field theoretic RG approach is to eliminate all singular 
ultraviolet cutoff dependences, i.e. singularities of the theory as (A, ^o) — > oo, by introduc- 
ing renormalization constants Zj, {i = 1,2, ■ ■ ■). RG flow equations are then obtained by 
examining how the Zi depend on the cutoff.lli In principle this approach is equivalent to 
the Wilsonian RG, which examines how the theory changes when the ultraviolet cutoff is 
changed from, e.g., A to A/6 with b a RG rescaling factor.li3 However, only a limited amount 
of work has been done on the formal relationship between these two formulations of the 

RG.ii 

A point of fundamental importance in any RG approach to a fleld theory is that one 
has to determine how many independent scaling operators there are. In the fleld theoretic 
RG approach, one needs to know how many renormalization constants are needed to make 
the theory flnite. The fleld theory deflned by Eqs. ( |2.lD -( P?BD has the form of a nonlinear 
sigma-model with perturbing operators. Let us flrst consider the pure sigma-model part, i.e. 
the first term on the r.h.s. of Eq. (|2.1|) with the constraints given by Eqs. (|2.3| ). This term is 
invariant under the symplectic symmetry group Sp(8nN), with N the number of replicas and 
2n the number of Matsubara frequencies. This model is well known to be renormalizable 
with two renormalization constants, one for the coupling constant G (the disorder) and one 
for the renormalization of the Q- fieldS The second term on the r.h.s. of Eq. (|2.1| ) breaks 
the symmetry to 

Sp{AN) ® Sp{4N) ® ■ • • ® Sp{4N) . 
-• , ■> 

2n times 

This term does not require any additional renormalization constants because the coupling 
constant H in Eq. ( |2.1| ) just multiplies the basic Q-field. and therefore the renormalization 
of H is determined by the field renormalization constant .011 This model represents the non- 
interacting localization problem.0 From a physical point of view the model has a rather 
restrictive property: The only interaction taken into account is the elastic electron-impurity 
scattering, and consequently the different Matsubara frequencies in Eq. ( pTl| ) are decoupled. 
Effectively, n is held fixed, and this is crucial for the simple renormalization properties of 
the noninteracting model. The situation remains relatively simple if further terms are added 
which respect the non-mixing of the frequencies. In this case, one needs one additional renor- 
malization constant for each operator that represents a different irreducible representation 
of Sp(8nN).§ 

The situation changes fundamentally with the addition of the last term in Eq. ( p.ll) . 
Physically, this term describes the electron-electron interaction and hence the exchange of 
energy between electrons. Technically, this leads to a coupling between the Matsubara fre- 
quencies, and an examination of the perturbation theory shows that this term introduces 
new infrared and ultraviolet singularities as n — > oo. In the absence of interactions, sin- 
gularities arise only from momentum integrations and the symmetry group is fixed to be 
Sp(8nN). With interactions there are singularities due to both momentum and frequency 
integrations, and the symmetry properties of the model change continuously during the RG 
procedure. As a consequence, the results mentioned above concerning the renormalizability 
of nonlinear sigma-models with perturbing operators, which apply to models with a fixed 
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symmetry, are inapplicable. No general results concerning the number of renormalization 
constants needed are available, and it is unclear how to renormalize the model given by the 
full Eq. ( |2.1| ). The symmetry arguments quoted above can provide only a lower bound on 
the number of Z's needed to renormalize the theory, and it is not known whether the model 
is renormalizable with a finite number of renormalization constants. 

All RG treatments of the field theory defined by Eqs. ( p.l|) -( ^73D so far have ignored this 
general renormalizability problem. They have assumed, explicitly or implicitly, that the full 
model is still renormalizable with one extra renormalization constant for each interaction 
coupling constant which is added. In addition, H acquires a renormalization constant of 
its own once interactions are present. In the absence of Cooperons, i.e. in a theory which 
contains only Kg and Kt, there is empirical evidence based on perturbation theory for this 
assumption being correct .011 In the presence of Cooperons things are more complicated as 
we will see. For pedagogical reasons, and to make contact with previous work, we neverthe- 
less proceed for a while using this assumption. It is most convenient to use a normalization 
point RG. The renormalized disorder, frequency coupling, interaction constants, and q-fields 
are denoted by g, h, fc^, fct, kc, and qR, respectively. They are defined by, 

G = ii'^Zgg , (3.1a) 
H = Zhh , (3.1b) 

Ks,t,c = Zs^t,cks,t,c , (3.1c) 

q = Z'/^qn , (3.1d) 

where n is an arbitrary momentum scale that is introduced in Eq. (|3.1a| ) to make g dimen- 
sionless. The renormalized N-point vertex functions are related to the bare ones by, 

rif )(p, Qn, 9, h, kt, k,- /i. A) = ^{N)^^^ jj^ _ (3 2) 

The theory is called renormalizable with these renormalization constants if all of the F^'' 
are finite as A — oo for fixed renormalized coupling constants. The Z's in Eqs. ( |3.1j ) are 
not unique. We fix them by the following normalization conditions for the two-point vertex 
functions, 

1 (2) f^^ 

rl^H )"jA"nin2(P)lp=0, w„j-w„2=M^M, a;„-^+a;„2=0 = —Sd/{'^T^) , (3.3a) 



j |'-p(2)\a/3,a/3 / xi 



(3.3b) 



1 



Ol-*- R )nin2,nin2\P ~ ^) Q \^ R Jnin2,nin2\P ~ 



= 27rTk^ 



(3.3c) 
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R )nin2,nin2\P ~ ^) 1 R Jnin2,ni?i2 ~ ^)\a^f3 



/p(2)NQ:/3,a/3 



-2nTkr 



(3.3d) 



These conditions determine the renormahzation constants Zg^h^s,t,c- The wavefunction or 
field renormahzation constant Z we fix by a normahzation condition for the one-point vertex 
function T^^\ We require, 



0/p(l)^Qa| _ 1 

R Jnn\Q„=fj.D/gh — 



(3.3e) 



Note that the normahzation conditions given by Eqs. ( p.3| ) are the usual ones: At scale /x 
the renormalized vertex functions are taken to have their tree level structure. 

From Eqs. (^), ( pJ[5| ), Q, and the Z's can be determined. We obtain, 



Z=(l/r«(fi„ = /x^M) 



(3.4a) 



9 (2vr 



+0(^^ 



(3.4b) 



(3.4c) 



■' s,t,c 



^s,t,c 



+0{g') , 



(3.4d) 



The fastest way to derive the RG fiow equations is to switch from our cutoff regularized 
theory to a minimal subtraction scheme. We can do so by putting e < in Eqs. 
letting A ^ oo, and then analytically continuing to e > 0. We find. 



4e 



— + 3lt + T 

e 



(3.5a) 



= 1 + £{5 - 3(1 + lht)k -l] + Oig 



^/^ = 1 + f -i + 37t + r 

8e L 



Zg — Zh 



Zt = l + - 
* 2e 



(3.5b) 

(3.5c) 
(3.5d) 
(3.5e) 
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1+f 

4e 



— + 3k + f 
e 



+ ^(l + 37j/7c 



(3.5f) 



with It = ln(l + 7t), 7f_c = kt,c/h. In giving Eqs. (|3.5|) , and below, we have for simphcity 
neglected terms of order gT"' with n > 2. The omission of these terms does not affect our 
conclusions. 

The one-loop RG flow equations follow from Eqs. ( |3.1| ) and ( |3.5| ) in the usual way.El With 
b ~ fi~^ and x = In b one obtains, 



dx 



-eg + ^{5-3{l + lht)lt-'-} + 0{g 



r- 



(3.6a) 



f = f [37*-l + f] + 0(/) , (3.6b) 

^ = |(l+7t)' + |(l + 7*)(l + 27*)f + 0(/) , (3.6c) 

^ = ?(37t + l) + ^{-- + 6/, + l-37. + f 
dx 4 8 e 

+--} + 0(/) . (3.6d) 

In Eqs. (p.5|), ( |3.6| ) T is the Cooper propagator at scale yu: 

f = f(/i) = f(^, n, uj') \^=^,=^D/2gh=n/2 ■ (3-7a) 
To lowest order in the disorder one has, 

fifi) = -^ + 0{g) . (3.7b) 
1 + 7cX 

In giving Eqs. (|3.5| ), (|3.6| ) we have specialized to the case of a Coulomb interaction between 
the electrons. In this case a compressibility sum rule enforces 7^ = kg/h = — 1,^ and the terms 
ln(H-7s) that appear in doing the integrals become — 2/e in Eqs. (|3.6|) . The presence of these 
terms in the RG flow equation for 7c reflects the well-known (In)^ singularity that exists in 
the disorder expansion of the single-particle DOS in D=2.i It has sometimes been claimed 
that these DOS effects are absent in all flow equations for the interactions constants.ii 
We stress that this statement depends on what quantity exactly one tries to derive a flow 
equation for. We find that it is true for 7^ and •jt, but not for the Cooperon amplitude 
7c. It is also important to note that the presence of these 1/e-terms in the flow equations 
per se does not create any problems. They do appear, e.g., in the renormalization of the 
single-particle DOS, for which a careful application of the RG^ leads to results that are 
consistent both internally and with those obtained by other methods.0 



15 



The Eqs. ( p.6|) are valid for the universahty class G. For the spin-orbit class SO the 
analogous results are, 

^ = -€<, + - f] + 0(/) , (3.8a) 

*._^|l_fl + 0(.^) , (3.8b) 

Note that the Eqs. (|3.8|) for the universality class SO have a fixed point describing a 
while the Eqs. (|3.6| ) for the generic universality class do not. The description of the MIT 
in the class G is more complicated. l3 For our present purposes we do not need a detailed 
description of the MIT and the Eqs. (|3.6|) are sufficient. 



B. A Scaling equation for T 



The RG flow equations given by Eqs. (|3.6| ) are not closed because they contain the 
Cooper propagator F. Ref. just used Eq. (|3.7bD in Eqs. (WM- In general this is not 
satisfactory, since for a consistent one-loop RG description one needs dV/dx to one- loop 



order. In Ref. |10| we argued that this was justified near the MIT since we found that 7c 
approaches a finite FP value at the transition. Here we try to improve on this point. In 
principle, F can be determined by solving the integral equation given by Eqs. (lOOp together 
with Eqs. ( p.6| ) and ( p.7| ). This structure is very different from the one usually encountered 
in RG approaches where the flow equations are an autonomous set of coupled differential 
equations. In this subsection we discuss attempts to reduce the integral equation for F to 
a single differential equation. We note that for generic integral equations this can not be 
done. What has effectively been assumed in the previous literatur^'i is that the reduction 
is possible in this case. As we will see, this reduction leads to severe structural problems 
which went unnoticed in the previous treatments quoted since the techniques used were not 
sensitive to them. In the next section we will therefore turn to a different method, which 
determines the behavior of F directly from Eqs. ( p. 201 ). This will qualitatively yield the 
same result as inserting Eq. ( 3.7b[ ) into Eqs. ( |3.(j| ). Here we investigate possible avenues for 



deriving a flow equation for F strictly in order to make contact with previous work. 

Within the normalization point RG approach we can formally obtain a differential equa- 
tion for F by renormalizing the Cooper propagator rather than the Cooper vertex function. 
We impose a normalization condition, 

f(/i)=7c , (3.9a) 

and define a renormalization constant Zc by 

7f = ^c7c . (3.9b) 
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Notice that this approach assumes that F is a simple scahng quantity. To zeroth order in 
the disorder Eqs. ( p. 71) and (|3.9| ) give, 

Z,= (l-7,x)-i , (3.10a) 



. (3.10b, 

A few remarks are in order in the context of Eq. ( p.lOb| ): (1) It has the standard form of 
a RG flow equation for a marginal operator. (2) In this approach the RG is used to obtain 
the BCS logarithm. This is in contrast to, e.g., Eq. (|2.11b| ) where we obtained the BCS 
logarithm at Gaussian order by inverting the vertex function. (3) In the present subsection 
the physical meaning of 7c is different from the rest of this paper. Due to Eq. ( p.9a| ) it plays 



the role of the Cooper propagator rather than that of the Cooper interaction amplitude 
that appeared in the previous subsection. (4) A crucial question is what the structure of 
the higher loop corrections to Eq. (|3.10b| ) will be. In a particular approximation that has 



been made in the literature and which we will discuss below, 7c flows to a nonzero fixed 
point value, 7c — >■ 7*- In this case Eq. ( p.lOa| ) naively implies that Zc has a Cooper-type 
singularity at a finite scale x = In 6 = I/7* which does not correspond to any physical phase 
transition. Of course, this conclusion is in general not necessarily correct, since terms of 
higher order in the disorder could change the behavior of Z^. Nevertheless, we will see that 
the appearance of such an unphysical singularity represents a serious problem. We will 
discuss this point in connection with Eq. (|3.14|) below, and in the Conclusion. 

Within this approach the one-loop RG fiow equation can be obtained by using Eqs. 
( p.l5| ) and (|2.19a|) in Eq. ( p.20a| ), and iterating to first order in the disorder. From Eqs. 



( ^.1| ) and (|3.9| ) one can then obtain a fiow equation for 7c. The resulting equation is quite 
complicated. Since we have come to the conclusion that this is not a viable approach we will 
not give a complete discussion, but rather illustrate only a few points. First, let us make 
contact with previous workEH by retaining only the first two terms in Eq. ( p.l5e|) , neglecting 



the corrections to H, putting p = 0, and working in D=2. In this approximation, 7 in Eq. 
( p.20a| ) is given to leading logarithmic accuracy by, 

7(a;,l^,^0 = 7f -ailn( ^+^^ + ^ ) , (3.11a) 



with, 



^ _ /(^/16)(-7s + 37t) for class G ^'i^^h^) 
^ \{g/16) for class SO ' ^ ' 



for the generic and spin-orbit universality classes, respectively. Using Eq. ( p.lla[ ) in the 
procedure discussed above, we obtain the following one-loop fiow equation for 7c, 



^ = -72[i _a^c] -4ai7cln2 + 2ai + 0(/,x2e-2^) , (3.12a) 



with. 
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|„f£±^-)_f! . (3.12b) 



z+1/2 \z + l/2J 12 



Two things should be noted: (1) The coefficients in Eq. ( p.l2a| ) are universal and, conse- 
quently, this equation has the form of a standard RG flow equation. (2) At the MIT in, 
e.g., the SO universality class, oi, which is of 0(g), is a constant of order g* = 0(e). If Eq. 
( ^.12a| ) is valid, then this implies that 7c goes to a flxed point value atthe MIT which is of 
order e^^^. In a strict e-expansion the terms of 0{g'~fc) and higher on the r.h.s of Eq. (p. 12a ) 



can be neglected, and to leading order the equation takes the form, 

--fj + 2ai . (3.12a') 



d^c .2 



dx 

This result is identical to those in Refs. |T] and |^, except for the prefactor of the ai, which 
is 4 in Refs. |l|,p. This difference is due to the fact that for a Coulomb interaction, which 
was considered in these references, an additional term appears in Eq. ( |2.15e| ). This leads to 
a more complicated kernel than the one in Eq. (|3.11a|) , and to an additional factor of 2 in 



Eq. ( p.l2a'D . Since this difference is not relevant for our purposes we restrict ourselves to 
the simpler kernel for the short-range case, Eq. ( |3.11a| ). 

In the approximation discussed above one flnds 7c const, at the MIT. As noted below 
Eqs. ( p.lO[ ), this result creates some problems. To illustrate this point we add to Eq. ( |3.11a{ ) 



via Eq. (|2.19a| ), the Cooper propagator contribution to -ffnin2 and Hn-^m, i-e. the integral 
over II in Eq. ( p^.l5b| ). The result is a kernel in Eq. ( |2.20aD which is given by. 



7(o.,l],^0 =7^°^ -ailn( "+^^+^ ) +a27f (ln[l-^7l°)ln(.;M)] 

+ ln[l-i7Wln(o.7fio)]) , (3.13a) 

with, 

aa = g/lQ . (3.13b) 
The last two terms in Eq. (p.l3aD occur because of the dependence of 7 on the Cooper 



propagator. In the last term in Eq. ( |3.13a| ) we have neglected a dependence on VL which is 



irrelevant for our purposes. Using Eq. p.l3a|) in the RG procedure yields the flow equation, 



^ = - aic - 2a,] - 4ai7cln2 + 2a, - a,^l ^^^^ ln[l - | ln(2^)] 



+0(/) . (3.14) 



The crucial point is that as x = In 6 ^ 00 the last term in Eq. ( |3.14|) does not exist because 



of a singularity at 6 ~ ea;p(l/7c) ~ exp{l/ e^^'^). Notice that this breakdown of the RG flow 
equation for the Cooper propagator is nonperturbative in nature. Previous attempts to de- 
rive a RG flow equation for the Cooper propagatoii'i were based on a frequency-momentum 
shell RG approach. This method is by necessity restricted to low order perturbative ex- 
pansions in both g and 7c and cannot be used to discuss the singularity in Eq. ( p. 141) 



Furthermore, if one expands the last term in Eq. (|3.14| ) in powers of 7c then a non-Borel 
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summable divergent series is obtained. The structure of this singularity resembles what is 
known as the renormalon problem in quantum field theory.0 We also note that this singu- 
larity is obviously related to the one discussed below Eq. ( p.lObp . 

In the following section we will argue that the Cooper propagator actually consists of 
multiple scaling parts, and that any theoretical approach that does not acknowledge this 
feature is invalid. It is unclear whether or not the singularity discussed in connection with 
Eq. ( |3.14D is related to this problem, which in turn is related to the question of whether or 
not the theory is renormalizable. Usually it is necessary to go to higher than one-loop order 
in order to verify the presence of multiple scaling parts. In the present case, however, the 
fact that a renormalization of the propagatorSS and the vertex function,^ respectively, led 
to different results is already an indication of their presence. 



IV. THE ELIASHBERG EQUATION, AND LOGARITHMIC CORRECTIONS TO 

SCALING 

In the first part of this section we discuss some general features of the Eliashberg equation 
with a repulsive kernel given by Eq. ( p.20a| ). We then give theoretical arguments in favor of 
the existence of logarithmic corrections to scaling at the MIT, and discuss their experimental 
relevance. 



A. The Eliashberg Equation 



In order to complete the RG description started in Sec. |111 A| , the scaling properties of the 
Cooper propagator given by Eqs. ( p.l9|) , (|2.20|) , and (p.7|) need to be determined. In Sec.|TTB 



we showed why previous attempts to reduce the integral equation for F to a differential 
equation are not satisfactory. Here we pursue a different approach. We acknowledge that 
one has to actually solve the integral equation in order to obtain information about F. Since 
this is very difficult to do in general, we classify possible solutions for different behaviors of 
the kernel 7(w, Q, uj'). 

To simplify our considerations we work at zero momentum. We further specialize to a 
model with a separable kernel. After drawing some conclusions for this special case we will 
argue that these are actually generic. The main points can be illustrated using a kernel that 
is a sum of two separable parts. 



^{00, n, 00") = Muj, n) Mlo", n) + Uu, n) Mu", n) . 

Note that Eq. (eq:4.1a) satisfies the symmetry requirement, 

7(c(j, Q, uo") = '~f{uj", uj) 



(4.1a) 



(4.1b) 



1712^13714 



K^n^m^ which in term 



Eq. ( |4.1b| ) is a consequence of the symmetry property K'^_^ 
follows from Eq. ( p.2c|) . Inserting Eq. ( |4.1a| ) into Eq. (p.20a|) leads to a separable integral 
equation that can be easily solved. We obtain. 
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t{uJ,n,Uj") = [(1 + Jl)(l + J2) - J; 



21-1 



X 



+f2{u,Q)f2{uj",Q)il + J,) 



with. 



r 

Jl,2 = / 

Jo 



2uj' + n 



and, 



^ Jo ^ 2u;' + n 



(4.2a) 



(4.2b) 



(4.2c) 



In the previous hterature it has been suggested that f (/i) given by Eq. ( p.7a| ) goes to a 
finite fixed point value at the MIT (cf. the discussion in the previous subsection), and that 
the approach to criticality is characterized by a conventional power law. To see how this 
type of behavior can be realized from Eqs. ( |4.1|) , (|4.2| ), we specialize to criticality, and put 
/2 = 0. If fi diverges like, e.g.. 



then the fixed point (FP) value of F is. 



, (a > 0) 



Sag— 2a 



dx 



(a; + l/2)(a;+ 1)2" 



-1 



(4.3) 



(4.4a) 



and near the FP F satisfies the fiow equation. 



dT 



d In /i^ 



t^2a/r -2af 



(4.4b) 



that is, F-F* ~ Similarly, if /2 = and /i vanishes at the MIT (Eq. (HJ) with a < 0), 
then F also vanishes and satisfies a fiow equation with universal coefficients. The marginal, 
logarithmic approach to zero occurs if /2 vanishes and /i approaches a constant at the MIT. 

If we assume that these asymptotic results are not tied to the separable kernel, but are 
generic properties of the general Eliashberg equation, then we have the following situation: 
If 7 diverges (vanishes) at the MIT, then F has a finite (zero) FP value, and if 7 goes 
to a constant then F vanishes logarithmically slowly. In all of these cases F satisfies an 
autonomous differential equation with universal coefficients. 

However, from a more general point of view, going beyond the asymptotic behavior, one 
expects a more complex result. For instance, even if 7 — 7* at the MIT one expects a 
correction that vanishes either as a power law or as a logarithm. In fact, the Eq. ( |3.6d[ ) 
predicts this kind of behavior. In our model calculation this happens if /2 7^ 0. With Eqs. 
( [4. 21) it is easy to see that F(/i) does not satisfy an autonomous differential equation if /2 7^ 0. 
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Of course, this just reflects the obvious fact that in general an integral equation cannot be 
reduced to a single differential equation. 

We further note that even if 7 diverges at the MIT, then one generically still expects a 
finite subleading contribution. Using this in either Eq. ([4.2a]) or in the actual Eliashberg 
equation, Eq. ( p.20a| ), one finds that (1) r(/i) does not satisfy a single differential equation, 
and (2) F approaches a finite FP value F*, but logarithmically slowly so. 

We conclude that for both the case where 7 diverges and where it approaches a constant 
at the MIT one expects a logarithmically slow approach to the FP value of F. The only other 
possibility is that 7 vanishes as a power law at the MIT. For this case F also vanishes as a 
power law, and in a scaling description (cf. below) F is a conventional irrelevant variable. 
While at present we cannot exclude this scenario, we consider it unlikely because of the first 
term on the r.h.s. of Eq. ( p.6d| ) (or the second term on the r.h.s of Eq. ( p.l5e| )), which tends 
to drive 7c and hence 7 towards larger values. 

B. Logarithmic Corrections to Scaling 

In the previous subsection we have argued that in general one expects F to approach 
its fixed point value logarithmically slowly. This result is consistent with Eq. ( |3.6d|) which 
gives 7c — i> 7c at the MIT, which in turn implies that F vanishes logarithmically slowly at 
the MIT. Because F couples to all physical quantities, cf. Eqs. (|3.6|), we conclude that in 
all universality classes where Cooperons are present, logarithmic corrections to scaling will 
appear. Note that this conclusion is independent of the spatial dimensionality and depends 
only on whether or not the kernel 7 has a constant contribution at the MIT. We also note 
that this is a zero-temperature, quantum mechanical effect that might be relevant for other 
quantum phase transitions. 

For a specific example of an observable quantity, let us consider the electrical conductivity 
cr, which is related to the disorder by, 

a = 8SDb-'/i2n)^7Tg{b) . (4.5) 

Note that in giving Eq. ([4.5| ) we have used units such that e^/^ is unity, and we have 
ignored the possibility of charge renormalization. For a discussion of the latter point we 
refer the reader elsewhere.E^ With t the dimensionless distance from the critical point at 
zero temperature, and 5r = F — F*, the conductivity satisfies the scaling equation, 

a{t, T) = J^[b^/^t, ¥T, 6r{b)] . (4.6) 

Here u is the correlation length exponent, z is the dynamical scaling exponent, and of the 
irrelevant variables in the scaling function JF we have kept only the one that decays most 
slowly at the MIT, i.e. (5F. 

At zero temperature we let b = t~'^, and assume that 0, x] is an analytic function of 
X since it is evaluated far from the MIT. We obtain, 

.(t^O,T^O)..otf + ^ + 

+ ...} , (4.7) 
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with s = v{D — 2). Here cxo is an unknown amplitude, and the are unknown expansion 
coefficients. Depending on what the subleading behavior of 5T{h) is, the Oj with i>2 could 
carry a very weak t-dependence (e.g. a2 ~ Inlnt). 

An interesting consequence of the logarithmically marginal operator F is that the dy- 
namical scaling exponent in Eq. (|4.6| ) is ill-defined. To see this we use that z is normally 
defined by, 

z = d + i^ , (4.8) 



with K the exponent that determines the anomalous dimension of h. The Eq. ( |3.6b| ) suggests 
that in general the linearized RG flow equation for h will have the form, 

d\\ih ^ ~ , ^ ^, 

_ — = K + a5T , (4.9) 

ax 

with a a universal constant and k = nii 5T vanishes as a power law. However, if 5T vanishes 
logarithmically at the MIT, then, 

h{h) = 6"(ln 6)" , (4.10) 

i.e., h does not scale as b'^. This in turn implies that Eq. ( |4.(j| ) should be replaced by, 

a{t,T) =b-'J^[b^/''t,b^h{b)T,6r{b)] . (4.11) 

At t = and as T ^ 0, the Eqs. (|4lO| ) and (gJT^) give. 



a{t = 0,T ^ 0) ^ (roT'/'[\n{l/T)Y''/' 

We conclude that for a(t = 0,T) the asymptotic scaling is in general determined by loga- 
rithms. 



V. EXPERIMENTAL RELEVANCE 

A. Experiments in Zero Magnetic Field 

It is well known that for the experimental determination of critical exponents, and for 
the comparison of theoretical and experimental values for these quantities, one must take 
into account corrections to scaling. This is so mainly because the asymptotic critical region, 
where corrections are negligible, is too small to be experimentally accessible. Furthermore, 
a reliable determination of the critical exponents usually requires experimental data that 



cover many decades of the control parameter.E3 For conventional phase transitions, where 
the control parameter is the temperature which is relatively easy to control, these conditions 
can be met. In the case of the MIT the situation is much less favorable. The main reason 
is that changing the control parameter, i.e. the impurity concentration, usually requires the 
preparation of a new sample. The only known way to avoid this problem is the stress-tuning 
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technique of Ref. Also, since the transition occurs at T = 0, measurements at very low 
temperatures and a careful extrapolation to T = are required. The application of these 
techniques to Si:P has led to the most accurate determination of the critical behavior of a 
MIT system to date.i@ Still, by the standards of critical phenomena experiments the data 
taken are relatively far from the critical point, with t > 10~^, and corrections to scaling 
have not been considered in the data analysis. This means that the measured exponent for 
the conductivity cannot be identified with the asymptotic critical exponent s. Rather, the 
measured value must be taken to represent some effective exponent Se//, which is different 
from s due to corrections to scaling. This observation offers an interesting possibility to 
explain the discrepancy between the measured value 0.51 for s, and the theoretical bound 
s > 2/3:i If the measured value represents an effective exponent, then the latter is not 
subject to the theoretical constraint for s. According to this hypothesis, the discrepancy 
between the measured value Sg// = 0.51 and the theoretical lower bound s > 0.67 must then 
be due to the corrections to scaling. With the usual, power-law corrections such a large 
discrepancy would be hard to explain. In the present case, however, where the corrections 
are logarithmic, it turns out that they provide a viable explanation for the observations. 

We can use Eq. ([4.7| ) to reconcile experiments!^ near the MIT which seem to give 
s < 2/3 with the theoretical boundi s > 2/3. In Fig. ^ we show experimental data, 
extrapolated to zero temperature, for the conductivity of Si:P. The data points were chosen 
as follows. For small t, roundoff due to sample inhomogeneities sets a limit at t ^ lO^^.iil 
At large t, at some point one leaves the region where scaling, even with corrections taken 
into account, is valid. We chose to include points up to t = 10~^. Obviously, for t ^ 1 the 
concept of corrections to scaling loses its meaning, and the expansion, Eq. ( [4.7| ), in powers 
of 1/lnt breaks down. We assumed a standard deviation of a quarter of the symbol size for 
all points except for the one at the smallest t, where we assumed it to be three times as 
large. In order to improve the statistics with many correction terms taken into account, we 
augmented the thirteen data points for 10"'^ < t < 10~^ by another twelve points obtained 
by linear interpolation between neighboring points. If all are set equal to zero then a 
best fit to the data yields s = 0.51 ± 0.05 < 2/3.i We now assume, somewhat arbitrarily, 
s = 0.7, and use Eq. ( [4.7[ ) with the ^ 0. We have used a standard fitting routine with 
singular value decomposition^ to optimize the values of the Oj. The dotted, dashed, and 
full lines, respectively, in Fig. |1] represent the best fits obtained with one, two, and three 
logarithmic correction terms. These fits are of significantly higher quality than a straight 
line fit optimizing s. More than three correction terms did not lead to further improvements 
in the fit quality. While the value of s was chosen arbitrarily, this demonstrates that this 
experient is certainly consistent with a lower bound s > 2/3 once the logarithmic corrections 
to scaling are taken into account. We have also tried to optimize s by choosing different 
values for s and comparing the quality of the resulting fits with a fixed number of coefficients 
taken into account (letting s float together with the proved unstable). The result was a 
very shallow minimum in the fit quality around s=0.7. Relatively large (±0.15) fluctuations 
in the best value of s were observed when large t data points were successively eliminated. 
From our fitting procedure for this experiment we estimate s = 0.70tom, 

where the lower 

bound is set by the theoretical bound rather than by the fit. 
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B. Experiments in a Magnetic Field 



While this success in fitting the Si:P data is encouraging, one might object that the model 
invoked contains an infinite set of unknown parameters, namely the in Eq. (eq:4.7), and 
that therefore the quality of the fit obtained is of no significance. Also, since the value of 
the asymptotic exponent in 3-D is not known, no quantitative statements can be made. 
It is therefore important to see if the model can predict any qualitative features that are 
independent of the actual value of s, and whether or not such features are observed. 

There are obvious features predicted by our model which follow from the qualitative mag- 
netic field dependence of the Cooperon. The first one is a strong magnetic field dependence 
of the effective exponent Sg//. If the logarithmic corrections to scaling are caused by the 
Cooper propagator, which acquires a mass in a magnetic field, then any nonzero field must 
act to destroy the logarithms. It thus follows from the model that Si:P, or any material, in a 
magnetic field must show a value of Se// that is larger than 2/3. Any observation of an Sg// 
smaller than 2/3 in a magnetic field would rule out the logarithmic corrections to scaling as 
the source of the anomalously small value of Se//. 

The second feature concerns the temperature dependence of the conductivity. It is well 
known that those materials which show an Se// < also show a change of sign of the 
temperature derivative da/dT of the conductivity as the transition is approached. Bill Within 
the RG description of the MIT this phenomenon can be explained by a change of sign of 
the 5f^-term in the g-fiow equation, Eq. ( |3.6a{ ) or ( |3.8aD . Since it is observed in Si:B, for 
which Eq. ( p.8a| ) is the relevant fiow equation, as well as in Si:P, the change of sign cannot 
be associated with the scaling behavior of the triplet interaction constant, but must be due 
to r. Since a magnetic field suppresses the Cooperon channel, this feature should therefore 
disappear in a sufficiently strong magnetic field. 

Let us make these predictions somewhat more quantitative. The relevant magnetic field 
scale is given by the magnetic length Ib = tic/eB being equal to the correlation length 
^ ~ kpt'". Let us assume, as we did in connection with Fig. 0, that the boundary of the 
critical region is given by t ~ 0.1. With kp 4 x lO^cm"^ for Si:B or Si:P near the MIT, 
and with i/ ~ 1 we then obtain ^ IT. The model thus predicts that for magnetic fields 
exceeding about IT the observed effective conductivity exponent should be larger than 2/3, 
and the change of sign of da/dT observed in zero field should disappear. 

Both of these features have been observed in the recent experiments by Sarachik et a/. .Bill 
The Cooperon induced logarithmic corrections to scaling thus provide a consistent explana- 
tion for several observed features of doped Silicon which otherwise would be mysterious. 



VI. CONCLUSION 

In this paper we have reconsidered the disordered electron problem in the presence of 
Cooperons, previous treatments of which had led to confiicting results. In particular we 
have analyzed the technical differences between Refs. which renormalized the Cooper 
propagator, and Ref. [1^, which renormalized the Cooper vertex function. Let us summarize 
the results of this analysis. 

(1) Crucial differences exist between the electron-electron interaction in the particle- hole 
channel and the particle-particle or Cooper channel, respectively. In the particle-hole channel 
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the conservation laws for particle number and spin enforce an essentially scalar structure 
of the propagators and vertex functions. As a result, vertex functions are simply scalar 
inverses of propagators, and standard renormalization procedures applied to either object 
lead to identical results. In contrast, in the Cooper channel there is no conservation law 
which would lead to an analogous simplification, and vertex functions and propagators are 
related by complicated matrix inversion procedures, cf. Sec. |I|. This leads to fundamentally 
different structures of the singularities in the two objects, 

which poses a problem for the renormalization, cf. Sec. [Ill A . 



(2) An attempt to reduce the Cooper propagator renormalization to a single RG flow 
equation encounters severe structural problems which we have discussed in Sec. [Ill B| . At the 
Gaussian level, the renormalizaton constant needed displays a BCS-like singularity at a finite 
scale. This leads to a renormalized vertex function which is not finite, but shows the same 
singularity. If one ignores this problem and proceeds to derive a flow equation for the Cooper 
propagator, the singularity in the Gaussian renormalization constant produces imaginary 
terms in the flow equation at one-loop level which resemble the renormalon singularity 
known in quantum field theory. This problem is nonperturbative in nature with respect 
to the Cooper interaction constant 7c. This explains why it went unnoticed in Refs. 
which expanded in both 7^ and in the disorder g. It invalidates the simple fixed point 
structure found in these references, according to which the Cooper channel interaction is a 
conventional irrelevant operator. 

(3) Any renormalization scheme for the Cooper vertex function must acknowledge the 
fact that the object that appears in perturbation theory is the Cooper propagator F, not just 
the Cooper interaction amplitude 7c. In fact the propagator plays the role of an effective 
interaction amplitude, cf. Sec. [Ill A . Due to the inversion problem mentioned in point (1) 
above it is a much more complicated object than 7c. This point was missed in Ref. |TU|. The 
procedure there was effectively equivalent to inserting the zero-loop order result for F into 
the one-loop order flow equations for the other coupling constants. The RG flow equations 
of Ref. |10| are therefore not consistent in general. They are valid near the transition only if 
7c has a finite FP value. 

(4) A possible solution of these problems is to derive flow equations for all coupling 
constants in terms of F, and to determine the latter from an Eliashberg-type integral equation 
which deals with the inversion problem. This is the approach taken in Sec. |IV|. Solutions 
of the integral equation obtained for separable model kernels suggest that F is not a simple 
scaling variable and does not satisfy an autonomous differential equation with universal 
coefficients. Plausible assumptions about the kernel of the integral equation lead to the 
conclusion that logarithmic corrections to scaling, as predicted in Ref. [TO, should exist 
regardless of whether F approaches a nonzero fixed point value, as asserted in Refs. ||,|], or 
vanishes logarithmically at large scales, as assumed in Ref. 0. However, this conclusion 
could not be made mathematically precise since the full integral equation has not been 
solved. 

(5) Logarithmic corrections to scaling, as discussed in Sec. IV B| , can explain some other- 
wise mysterious properties of materials that belong to universality classes with Cooperons, 
cf. Sec. |V]. The discrepancy between the asymptotic critical exponent for the conductivity 
and the observable effective exponent in this case is large enough to provide an explanation 
for the 'exponent puzzle' in Si:P. This explanation is bolstered by two predictions which 
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have been verified by recent experiments:^ In the presence of a sufficiently strong magnetic 
field the observed conductivity exponent s must satisfy the lower bound s > 2/3, and the 
characteristic change in the temperature dependence of the conductivity observed in zero 
field as the transition is approached should disappear. 

While these results are encouraging, some serious problems remain, most notably the 
question of whether or not the theory is renormalizable. We recall (cf. the discussion in 
Sec. [Ill Aj) that the pure nonlinear sigma-model is renormalizable with two renormalization 



constants. The proof of this makes heavy use of the symmetry properties of the model.ill 
Since the interaction terms in the action break the symmetry of the sigma-model, the proof 



of Ref. ^ is no longer applicable for the interacting model and renormalizability has never 
been proven. It is conceivable, however, that the conservation laws in the particle-hole 
channel are still sufficient to ensure renormalizability for the model without Cooperons. This 
would be consistent with the considerable body of perturbative evidence of renormalizability 
for this case. For the case with Cooperons the absence of a conservation law and the 
unusual structure encountered in attempts to apply renormalization techniques makes one 
much more doubtful. It would be useful to further study this question. For instance, it 
would be interesting to see whether the conservation laws are indeed sufficient to ensure 
renormalizability. It would also be useful to have an alternative approach to the subject 
which avoids the generalized nonlinear sigma-model with its many technical problems. One 
could, for instance, try to work directly with the Grassmannian action as in recent work 
on interacting fermions without disorder.!! One could also imagine an order parameter 
description of the MIT, using the Q-field theory, but not the sigma-model. Work in these 
directions is underway and will be presented in future publications. 
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FIGURES 



FIG. 1. Zero-temperature conductivity of Si:P vs. electron concentration. The data 
points have been redrawn from Ref. ^. The dotted, dashed, and full lines, respectively, 
are best fits using one, two, and three corrections terms in Eq. ( |4.7D with an asymp- 
totic critical exponent s = 0.7. Best fit values for the coefficients in Eq. ( |4.7| ) are 
do = 54.87; 94.82; 132.16; ai = -1.84; -4.42; -6.61; 02 = 0; 6.21; 17.73; 03 = 0;0; -16.38 for the 
three curves, respectively. See the text for more details. 
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